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Abstract. Kinetically constrained spin models are schematic coarse-grained 
models for the glass transition which represent an efficient theoretical tool to study 
detailed spatio-temporal aspects of dynamic heterogeneity in supercooled liquids. 
Here, we study how spatially correlated dynamic domains evolve with time and 
compare our results to various experimental and numerical investigations. We find 
that strong and fragile models yield different results. In particular, the lifetime of 
dynamic heterogeneity remains constant and roughly equal to the alpha relaxation 
time in strong models, while it increases more rapidly in fragile models when the 
glass transition is approached. 



1. Introduction 

The viscosity of supercooled liquids increases extremely rapidly when the temperature 
is reduced towards the glass temperature. It is firmly established that this dramatic 
slowing down is spatially heterogeneous. Local relaxation is widely distributed in 
time — existence of broad stretched relaxations — , but also in space — existence of 
dynamic heterogeneity pp. The main physical aspect is that spatial fluctuations of 
local relaxations become increasingly spatially correlated when temperature decreases. 
Direct experimental investigations of the time and temperatures dependences of the 
relevant dynamic lengthscales at low temperature are however still missing. 

To study dynamic heterogeneity, statistical correlators which probe more than 
two points in space and time have to be considered ^[21 EH- For example, if one wants 
to study spatial correlations of the local dynamics one has to define a two-point, 
two-time correlator, 

C 2 , 2 (N - j\,t) = <P(0, t)P 3 (0, t)) - (Pi(0, t))(P,(0, t)>, (1) 

where notations are adapted to lattice spin models. In equation QJ, denote lattice 
sites, Pi(0,i) quantifies the dynamics at site i between times and t (autocorrelation 
or persistence functions), and brackets represent ensemble averages. The physical 
meaning of Q is clear: given a spontaneous fluctuation of the two-time dynamics at 
site i, is there a similar fluctuation at site j? The quantity Q has now been discussed 
both theoretically and numerically in some detail [3J0J|S], generically revealing the 
existence of a growing spatial range of dynamic correlations in supercooled liquids 
accompanying an increasingly sluggish dynamics. 
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Logically, the next question is then: given spatial structures of the local relaxation 
between times and t, what will this structure look like, say, between times t and 2i? 
In other words (i , how do dynamic heterogeneities evolve with time? This question is 
in fact simpler to address experimentally because no spatial resolution is needed and 
different experimental techniques can be devised: NMR ;T, solvation dynamics [5], 
optical |3 and dielectric ^H] hole-burning. In statistical terms, one wants to study a 
four-time correlation function of the general form 

C 4 (h,t w ,t 2 ) = {Pi(Q,ti)Pi{t 1 +U»ti+U l +t2)), (2) 

which correlates dynamics between times and ti and between times t\ + t w and 
ti+t w +t 2 . Again the physical content of @ is clear JT]: given a dynamic fluctuation 
at site i between in a certain time interval t%, how long does it take for this fluctuation 
to be washed out? This leads to the general concept of a lifetime, Tdh, for dynamic 
heterogeneity. While many investigations ^ [21 EI d3 El Q3| indicate that Tdh is 
in fact slaved to the alpha- relaxation time of the liquid, Tdh ~ T a , photobleaching 
experiments very close to the glass transition indicate that Tdh may become several 
orders of magnitude larger than r Q , although in a surprisingly abrupt manner |14j . 

In this paper we study the lifetime of dynamic heterogeneity in kinetically 
constrained spin models of supercooled liquids These models represent schematic 
coarse-grained models for the glass transition and provide a very efficient tool to 
study in detail many spatio-temporal aspects related to dynamic heterogeneity such 
as dynamic lengthscales ^B], scaling JJjj or decoupling phenomena ^JEIj- They 
are simple enough that analytical progress can be made and numerical simulations 
performed on a wide range of lengthscales and timescales, and yet rich enough that 
direct comparisons to both simulations and experiments can be made. 

2. Models 

Following previous works El El El I2U| . we focus on two specific spin facilitated 
models in one spatial dimension, namely the one-spin facilitated Fredrickson- Andersen 
(FA) model [21] and the East model [22] that respectively behave as strong and fragile 
systems |15j . These are probably the simplest models which incorporate the ideas 
that (i) mobility in supercooled liquids is both highly localized and sparse, as revealed 
by simulations |12) : (ii) a localized mobility very easily propagates to neighbouring 
regions, the dynamic facilitation concept. Detailed studies in spatial dimensions 
larger than one have shown that dimensionality does not play a relevant qualitative 
role ^lEDGH], and justify therefore the present one-dimensional studies. 

Both models are defined by the same non-interacting Hamiltonian, H — ^i 71 *! 
expressed in terms of a mobility variable, m — 1 when site i is mobile, m = otherwise. 
Dynamic facilitation is incorporated at the level of the dynamic rules through kinetic 
constraints. In the FA model, the site i can evolve with Boltzmann probability if at 
least one of its two neighbours is mobile, + n i+ i > 0. In the East model the site 
i can evolve only if its left neigbour is mobile, rii_i = 1. 

We have performed numerical simulations of both models using a continuous time 
Monte Carlo algorithm where all moves are accepted and the time is updated according 
to the corresponding statistical weight. Simulations have been performed 'only' over 
about 7 decades in time because extensive time averaging is required to accurately 
measure multi-time correlation functions such as equation |2|). 
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Figure 1. Left: Four-time correlation function C^,(r a , t w , £2) normalized by 
its t 2 = value for various t w = 0, 550, 2100, 5200, 20000, 49500, and 191000 
(from right to left) in the FA model at T = 0.3. These dashed curves converge 
at large t w to the bulk persistence function, p{t2), shown as a full line. Right: 
Corresponding logarithmic distributions of relaxation times. 



3. Results 

3.1. Dynamic filtering 

There are several parameters involved in the four-time correlator J5J that need to be 
appropriately chosen. Since dynamic heterogeneity is more pronounced for times close 
to r a it is sensible to first fix t\ — r a and to study the remaining t w and ti dependences. 
As a local dynamic correlator we first focus on the persistence function [23|, Pi(0, t) = 1 
if spin i has not flipped in the interval [0, t], Pi(Q, t) = otherwise. We also define the 
mean persistence, p(t) = (Pi(t)), from which we measure r Q via p(r a ) = e . 

In figure ^ (left) we show the £2 dependence of Ci{T a ,t Wl t2) for various i w at 
T = 0.3 in the FA model. We have normalized C4 by p(r a ) — e _1 , its value at 
ti = 0. By definition, this function describes the persistence function in the interval 
[tw + T a, tw + T a + t%\ of those sites which had not flipped in the interval [0, T a ], and 
were therefore slower than average. The first term in the correlator J2J) plays the 
role of a dynamic filter , selecting a sub-population of sites which have an average 
dynamics different from the bulk. From earlier works studying the spatial correlator 

it is known that those sites belong to compact clusters that represent the largest 
regions of space with no mobility defects at time 1171 1^3*] . 

Immediately after filtering one expects therefore those slow regions to remain 
slow, as indeed observed in figure for t w = 0. When t w increases, this selected 
population gradually forgets it was initially slow. When t w — > 00, bulk dynamics is 
recovered, 

Ci{t\, t w — > 00, i 2 ) , . , , 
; — \ >P{t2), (3) 

P(Ta) 

as demonstrated by the full line in figure ^ I n figure (right) we also show the 
(logarithmic) distribution of relaxation times corresponding to the functions shown in 
the left panel, a representation sometimes preferred in experimental works |J. Both 
quantities are of course fully equivalent |20]. It is clear from these figures that once a 
subset of sites has been dynamically selected the remaining relaxation is narrower than 
the bulk relaxation. In fact all persistence functions shown in figure^are well described 




Figure 2. Fs(k,t) and its 'heterogeneous', F s (k,t), and 'homogeneous', 
[F a (k, t/2)] 2 , limits. Left: East model for fixed k = n/5 and various temperatures. 
Right: East model at T = 0.4 and various wavevectors. 



by stretched exponentials. While (3 = 1/2 is observed for the bulk dynamics, one finds 
P ~ 0.83 at t w = 0. Accordingly, distribution of relaxation times progressively broaden 
when t w increases. These results are consistent with experimental observations. 

In the FA model their interpretation is straightforward. Stretching in this model 
follows from an exponential distribution of distance between mobility defects jlfij . 
Dynamic filtering implies that this domain distribution is cut-off at small distance. 
Narrower lengthscale distributions directly imply narrower timescale distributions. 

We have also investigated the effect of changing the 'filter efficiency' .?, which 
in our case implies changing the duration of the filtering interval, [0, £i]. While bulk 
distributions are found for ti/r a <C 1 (weak filtering) distributions shift to larger times 
and become very narrow when t\jr a increases. In the following we work at constant 
filtering, t\ — r a . 



3.2. 'Homogeneous' vs. 'heterogeneous' dynamics 

The ability to select a sub-ensemble of sites that are slower than average is sometimes 
taken as a definition of dynamic heterogeneity 0, although lengthscales play no role 
in this view. That FA and East model display spatially heterogeneous dynamics is 
well-known, and the results of the previous section are therefore natural. 

Another indicator of dynamic heterogeneity has been proposed [111 1251 125] based 
on the analysis of the four-time correlation (0). Consider the situation where t w = 0, 
and t-y = t2 = t/2. In that case, one studies a 'three-time' correlation 

F 3 (t) = (P i (0,t/2)P i (t/2,t)). (4) 

Two extreme behaviours can be expected for F 3 (t). (i) Dynamics in the intervals 
[0,i/2] and [t / 2 , £] are totally uncorrelated, and thus F 3 (t) sa [p(t/2)] 2 . (ii) Dynamics 
in the two intervals are strongly correlated, in the sense that those regions that survive 
filtering in [0,i/2] are also those dominating the relaxation in the full interval [0,i]. 
In that case, Fa(t) as pit). Scenarii (i) and (ii) have been termed 'homogeneous' 
and 'heterogeneous', respectively although again lengthscales play no role in the 
distinction. Clearly both estimates become equivalent when p(t) decays exponentially. 

Of course when studying the persistence function in the FA and East models, 
scenario (ii) strictly applies by definition, because Pi(0,t/2)Pi(t/2,t) = P,(0, i). In 
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real materials, smoother dynamic functions are studied, directly defined from the 
particles positions instead of a mobility field. Our strategy is therefore to couple 
probe particles to our mobility field, see Refs. ^1 ^5] for technical details. From 
probe molecule displacements, Sx(0, t) — x(t) — x(0), we define self-intermediate 
scattering functions, F s (k,t) — (cos[fc ■ Sx(0,t)]), and the analog of equation J3J|, 
F 3 (k,t) = (cos[k ■ dx{0,t/2)]cos[k ■ Sx(t/2,ty\). 

Our numerical results are presented in figure Clearly the time dependence of 
F3 closely follows the one of F s (k,t), in agreement with the 'heterogeneous' scenario 
described above. This is consistent with numerical results 

In the present approach, this result is a natural consequence of decoupling 
between structural relaxation and diffusion E3 At large wavevectors, k ~ n, 
corresponding to distances of the order of the lattice spacing, F s (k, t) is dominated by 
the time distribution of the first jump of the probe molecule in the interval [0,i], 
so that F s (k, t) w p(t) w F 3 (k,t). At large distance, k < k* , Fickian diffusion 
holds 19 , F s (k,t) = exp(-fc 2 L> s i), and there is no distinction between homogeneous 
and heterogeneous relaxation. At intermediate wavectors, 7r > k > k*, the long-time 
decay of F s (k, t) is again dominated by the persistence time distribution, because the 
timescale it takes a molecule to make 2-x/k steps is strongly dominated by the timescale 
to make the first step This is just the condition for the heterogeneous scenario 
to hold, in agreement with figure|21 The characteristic wavevector separating the two 
regimes, k*(T) = l/*jT a D 8 , decreases when temperature decreases, opening a larger 
heterogeneous window; k* also sets the upper limit of validity of Fickian diffusion in 
supercooled liquids [T3| . 

3. 3. Lifetime of dynamic heterogeneity 

After dynamic filtering it takes some time for filtered distributions to reequilibrate 
towards the bulk relaxation, cf figure ^ To extract the typical lifetime of dynamic 
heterogeneity, Tdh, we tried several procedures which all lead to similar results, 
based on how timescales (time decay of persistence functions or moments of the 
corresponding distributions) return to their equilibrium values. Following Refs. |13II18| 
we also measured the integrated difference between filtered and bulk dynamics, 
A(t w ) = dt 2 [C4(T a ,t w ,t 2 )/p(T a ) - p(t 2 )}. From figure [TJ we expect that A(t w ) 
goes to on a timescale Tdh- In practice, we define Tdh as A(Ydh)/A(0) = e _1 . In 
principle, Tdh depends on the filtering time, t\, and on temperature, T. 

We show in figure E3 results at various T but constant filter efficiency, t\ = r a (T), 
in East and FA models. While Tdh is set by r a in the FA model (the tiny deviation 
observed in figure [3] is due to finite T corrections which weaken when T gets lower), 
this is not true in the East model where Tdh systematically grows faster than r a at low 
T, as emphasized in the inset. Quantitatively a power law relationship, Tdh ~ r a + ^i 
with C w 0.06, is a good description of the data, although alternative fitting formula 
could probably be used. 

In the fragile case, Tdh can therefore be considered as an additional slow timescale 
characterizing the alpha-relaxation |14|. on top of t q and 1/D S |18| . The comparative 
study of FA and East models offers a possible physical interpretation. While both 
models display stretched relaxations, in the FA model stretching is constant, (3 = 1/2, 
while (3 increases linearly with T in the East model [201 - Therefore r a represents the 
first moment of a distribution that becomes wider and wider when T decreases. We 
attribute the small but systematic decoupling between Tdh and r a to this broadening. 
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Figure 3. Integrated difference between filtered and bulk dynamics at constant 
filtering time, t\ = T a , and various temperatures in the FA (left) and East 
(right) models. Times have been rescaled to collapse the data points at various 
temperatures and extract the lifetime of dynamic heterogeneity, r^h- In the FA 
model, is set by r a with tiny corrections that vanish at low T. In the 
East model, T$h grows faster than r a , as emphasized in the inset showing the 
systematic increase of r^/ra as T decreases. Mastercurves are fitted by stretched 
exponentials shown with full lines, (3 = 0.75 (left) and (3 = 0.8 (right). 



Unfortunately this decoupling does not quantitatively account for the results of 
photobleaching experiments which show that T^h/T a increases strongly close to T g |14j . 
In OTP, while r a changes by about 1 decade when T is changed from T g + 4 K to 
Tg + 1 K, the ratio Tdh/i~ a changes by 2 orders of magnitude, so that («2. This value is 
much too large to be accounted for by the above results. Presumably, also, (3 does not 
vary much on such a tiny temperature interval. Therefore, the present results cannot 
explain the experimental value C ~ 2 without invoking possible non-equilibrium effects 
due to the proximity of T g . However we were able to predict instead a smaller, but 
definitely non- vanishing decoupling between the lifetime of dynamic heterogeneity and 
the alpha-relaxation time which could be detected in dynamic filtering experiments 
performed on a sufficiently large temperature window in fragile glass-formers. 
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